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ABSTRACT 

It has been suggested that multiple scattering on circumstellar dust could explain the non-standard 
reddening observed in the line-of-sight to Type la supernovac. In this work we use Monte Carlo 
simulations to examine how the scattered light would affect the shape of optical lightcurves and 
spectral features. We find that the effects on the lightcurve widths, apparent time evolution of color 
excess and blending of spectral features originating at different photosphcric velocities should allow for 
tests of the circumstellar dust hypothesis on a case by case basis. Our simulations also show that for 
circumstellar shells with radii r = 10 16 — 10 19 cm, the lightcurve modifications are well described by 
the empirical Am 15 parameter and intrinsic color variations of order <Jbv = 0.05 — 0.1 arise naturally. 
For large shell radii an excess lightcurve tail is expected in B band, as observed in e.g. SN 2006X. 
Subject headings: dust, extinction — Interstellar medium, nebulae :distance scale — Cosmology : 
supernovae: general — Stars : supernovae individual SN 2006X - Stars 



1. INTRODUCTION 

A dramatic breakthrough in cosmology took place 
more than a decade ago when the accelerated ex- 
pansion of the universe was discovered through 
Type la su pernovae (SNe la) used as standardiz - 
able candles (|Riess et all 119981: iPerlmutter et all 119990 . 
SNe la remain among the best tools to investi- 
gate the energy content of the universe, as demon- 
strated by recent v e ry su c cessful surveys of high- 
z SNe (lAstier et al.1 120061 iWood-Vasev et all [20071: 



Riessetal.l l2007t iKowalski et all 120081: iKes sler et al I 
20091 : lAmanullah et al.ll2010D . However, the progress for 
precision cosmology is hampered by systematic uncer- 
tainties, notably a potential drift in the SN la luminosity 
and intrinsic colors with cosmic time and the impact of 
changing dust extin ction properties along the line of sight 
(jNordin et al.ll2008[ ). The standardization of SNe la gen- 
erally uses two parameters: the measured color, typically 
rest-frame B — V and often one more optical color, and 
an empiricall y defined lightcurve shape parameter, such 
as st retch, s flPerlmutteilll997t iGuvl 120051: iConley et al] 
20081). xi (IGuv et all I20 07}h Am^(B ) (IPhillipsI 11993 
Hamuv et all 119961) or A ([Riess et all 119961: Uha et al.l 



2007TT In this work we examine how multiple scatter- 
ing of supernova photons on the surrounding material 
could affect both the lightcurve shapes and colors. 

While the reddening of SNe la was originally thought 
to be entirely due to extinction by interstellar dust, there 
is an increasing body of evidence showing that the mea- 
sured color excesses to some SN la show a steeper wave- 
length dependence than what is observed for reddened 
stars in the Milky Way. Likewise, there is no generally 
accepted model that explains the lightcurve shape vari- 
ations observed in SN e la and its co rrelations with peak 
supernova brightness (lPhillipslll993l). 

Following up the work by I Wand ()2005[ ). [G^oban (|2008l ) 
(hereafter G08) showed that multiple scattering on cir- 
cumstellar dust could potentially help explaining the 
low values of the total-to-selective extinction param- 



eter, R v = A V /E(B - V) < 3, observed in the 
sight lines o f near-by Type la SN e (lElias-Rosa et all 
20061 [20081 IKrisciunas et all [20071 IWang et al.l T2008F 



Nobili fc Goobar!l2008J ). 



A Monte Carlo ray tracing technique was used in 
G08 to investigate the reddening effects in the pres- 
ence of circumstellar dust and an approximate redden- 
ing law for opt i cal an d near-IR wavelengths was derived. 
iFolatelli et al.l (|2010l) compared the parametrized model 
in G08 with high quality uBgriY JHK data from the 
Carnegie Supernova Project and found that it provided 
particularly good fits for the two most reddened SNe la 
in their sample, SN 2005A and SN 2006X. The success of 
these comparisons calls for further scrutiny of the model. 
In this work we investigate the impact of multiple scat- 
tering on the shapes of SNe la optical lightcurves. We 
explore how the extra wavelength dependent "random 
walk" of photons in the presence of scatterers affects the 
width of broad-band filter lightcurves. In particular, red- 
dened "slow decliners" may be compatible with the sce- 
nario in G08. 

Furthermore, also the strengths of spectroscopic fea- 
tures as a function of time are likely to be affected, since 
a varying path-length would lead to a superposition of 
several SN phases. Throughout the paper we point out 
how observations could put bounds on the size of a hy- 
pothetical circumstellar dust shell. 

2. CIRCUMSTELLAR DUST IN THE SINGLE 
DEGENERATE SCENARIO 

In the currently favored Type la scenario, a white 
dwarf accretes mass from a companion star until the 
Chandrasekhar limit is reached, at which point insta- 
bilities lead to an explosion. However, the model re- 
quires the accreting white dwarf to expel material ac- 
cumulating on its surface through fast stellar winds, 
v w ~ 10 2 — 10 3 km s _1 , to avoid the formation of a com- 
mon envelope between the white dwarf and the donor 
star which would prevent the explosion (lHachisu et aLI 
119961: lHachisu fc Katolll999l: lHachisu et allll999fl . Con- 
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densation of dust grains may take place in the ejected 
matter, possibly also in collisions with the accreting mat- 
ter. For dust grain growth in stellar winds, the density 
as a function of distance to the star surface, r, is generi- 

cally expected to scale as pdust oc , where, v w , is the 

stellar wind velocity and M w is the mass-loss per unit 
time, typically of order 10 -7 M Q yr" 1 in SNc la models. 

3. CIRCUMSTELLAR DUST DESTRUCTION BY SN la 
LUMINOSITY 

So far, we have assumed that the supernova explosion 
site is surrounded by dust. However, the radiation from 
the supernova itself will sublimate dust grains as these 
are heated up to T su b ~ 3000 K. Thus, assuming that 
the entire UV-optical luminosity of the supernova, is ab- 
sorbed by surrounding dust grains, we can compute the 
radius, r c , where dust grains, if present at the explosion 
time, would be depicted. From the radiation balance 
condition: 



Lbolif) Qabs 

l6iraT sub < Q IR > 



1/2 



(1) 



where Q a bs is the effective absorption efficiency and 
< Qtr > is the Planck- ay eraged emissivity at T = T su b- 
IWaxman fc Draine! (|2000|) studied a similar situation, 
but for dust surrounding gamma-ray bursts and found 
that for 2000 < T < 3000 K, < Q IR >sa 0.1 for a 
mixture of astronomical silicate and graphite gr ains with 
typical sizes of 0.1 f»m. iSollerman et al.l ([2004D have es- 
timated the peak bolometric luminosity of SN la to be 
L bo i = 1.43 • M Ni ■ 10 43 ergs/s, where the 56 Ni mass 
(in units of solar masses), M Ni , has been mea sured to 
be between 0.1 and 1.1 M Q (jLeibundgutl 120001 ) . Insert- 
ing this peak luminosity (E 1 ^1 — 5 eV) into Eq. (TTJ), 
and assuming that all the radiation is absorbed, we find 
fc?S10 16 cm (~ 0.003 pc), which we will use as a rough 
estimate of the region where evaporation of pre-existing 
dust would take place as a result of the SN heating. A 
similar estimate of t he region that ought to be depleted 
of dust was done by iPearce fe MavesTl) 19861) . Thus, to 
create a dust shell that would survive the radiation from 
the supernova explosion, and using the stellar wind ve- 
locities and mass-loss rates from Section [2j the wind 
must start about 10-100 years prior to the explosion, i.e. 
J>1O -6 M0 of mass must be expelled for any effect of cir- 
cumstellar dust to be noticeable. For SN 2005gj, a well- 
studied example of a supernova with sig nificant interac- 
tion w ith the circumstellar environment, lAldering et al.l 
(|2006|) find evidence for mass-loss about 8 years before 
the explosion. For significant opacity, t ~ 1, the re- 
quired dust mass ca n be estimated from the ab sorp- 
tion cross-sec tions in iWeingartner fc Draind (|2001[ ) and 
IDraind (|2003t) . a a /m dust ~ 10 4 cm 2 /g- Thus, for a thin 
dust shell at r ~ r c ~ 10 16 cm, this corresponds to 
mdust ~ 1O~ 4 M . 

We also note that since r c oc y/Mm, intrinsically 
brighter explosions are likely to have larger r c , i.e., thin- 
ner circumstellar shells. In the following sections we will 
see that, depending on the size of the outer shell, mul- 
tiple scattering could generate relations between bright- 
ness and lightcurve shapes and colors qualitatively sim- 
ilar to the observed ranges. We note, however, that in- 
teraction with dust in the CS material can only broaden 



the pristine lightcurve, i.e., it is not possible to produce 
a faster decline than what results from the original su- 
pernova explosion. 

4. MULTIPLE SCATTERING OF PHOTONS AROUND THE 
SUPERNOVA 

A natural effect of light scattering is that it adds 
flight time to th e photons, t hereby affecting the ob- 
served lightcurve (jPatatll2005t IPatat et al.ll2006D . Next, 
we start by generalizing the iGoobarl ( 20081 ) model: we 
consider a scenario where the exploding white dwarf 
(rwD ~ 10 9 cm) is at the center of a homogeneous dust 
shell with an outer radius, r (r 3> ^wd), and an inner 
radius r;, where rwD < t% <t. Note that this will not 
change any of the results in G08, where ?*j = rwD ~ 
was assumed, but it will affect the distribution of photon 
flight times, i.e., the subject of this paper. 

The technique in G08 is extended to investigate the 
time domain by calculating the total traveled path 
length, D, for each photon before it crosses a plane, P, 
that is perpendicular to the line of sight between the ob- 
server and the source, and is located at the the dist ance, 
r, from the source which is illustrated in Fig. IA1I The 
details for the c alculation of t he photon travel times are 
outlined in Sec. |Appcndix A[ The traveled photon path 
length, D, is given by Eq. I All 

The traveled photon path length is distributed in the 
interval r < D < oo, with D = r for photons that do not 
interact with the dust at all. The shape of the distribu- 
tion of D will depend on the radius, r, the thickness of 
the shell, the scattering and absorption properties of the 
dust, and the dust density. In this work we will quantify 
the latter in terms of the observed color excess, E(B—V), 
relative to the pristine source color. For the scattering 
and absorption properties we use models consisting of 
carbonaceous grains and amo rpho us silicate gra ins from 
IWeingartner fc Draind ([20(51 and IDraind ([200l . 

We begin by exploring the scattering distribution for 
the two extreme scenarios (a) = rwD ~ (i.e., a ho- 
mogeneous dust sphere as in G08), and (b) r» = 0.95 • r, 
i.e., a very thin dust shell. Clearly, the scenario (a) is 
already challenged by the conclusions in Section [3l nev- 
ertheless, it serves as a point of reference. The left panel 
of Fig.[T]shows the distribution of D/r from Monte Carlo 
simulations of photons with A = 5000 A. Most photons 
travel unscattered, and for the rest, there is a rough ex- 
ponential fall-off of delay times. The dust density was 
chosen to cause an observed color excess E(B — V) m 0.5. 
Furthermore, we are assuming that the dust grain prop- 
erties correspond to LMC-type dust, although the time 
distributions for MW-type dust are very similar. 

A local minimum can be noted for case (b) around 
D/r ~ 2. This can be understood by realizing that pho- 
tons from the source can only double the path length, 
D/r ~ 2, by having a total scattering angle of ~ ir/2 
for this particular setup, which will take it through more 
dust than for other angles. It is then likely to scatter 
again, which would result in either a shorter or a longer 
path. The significance of the minimum will decrease with 
the density of the dust and the inner radius of the shell . 



LIGHTCURVE PERTURBATION FROM MULTIPLE 
SCATTERING 
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Figure 1. Monte Carlo results of two extreme circum stellar dust scenarios with a dust density corresponding to an observed E(B — V) = 
0.5. The upper and lower panels show the results for n = 0.0 and n = 0.95 • r respectively. The left panels show the histograms of the 
total traveled photon paths D (sec Sec. | Appendix A] for definition) normalized with the radius, R, of the sphere, while the middle panels 
show the histogram of the total scattering angles, and the right panels show the number of scatters. The local minimum around D/r ~ 2 
for the lower left panel is related to the lack of photons in the lower middle panel around cos(Stot) ~ with n sca t = 1. 



Since the scattering cross-section, scattering angle and 
albedo are wavelength dependent, the lightcurve pertur- 
bations are expected to vary between broad-band filters. 

We use the time dependent s pectral energ y distribution 
(SED) SN la templates from lHsiao et all (|2007t ) (here- 
after H07) to build up observed broad-band lightcurves 
in the presence of circumstcllar dust. Thus, we make 
the simplifying assumption that the empirically derived 
spectral template is pristine, i.e., not already affected 
by dust. This is not a major limitation for this analy- 
sis, since we are here primarily concerned with relative 
quantities, such as time-delays or color excess rather than 
absolute times or fluxes. 

Li ke in G08, we only conside r d ust grain mo dels 
from lWeingartner fc Draind (|2001D and lDraind (|2003D for 
Milky- Way dust and LMC dust. In particular, it was 
found in G08 that LMC dust grains in the circumstcl- 
lar medium would lead to Ry ~ 1.65, while Milky- Way 
dust lead to Ry ~ 2.56. These two grain models ap- 
proximately bracket the empirically found values for Ry 
for a large fraction of SNe la. In both cases, the extinc- 
tion law deviated sig nificantly from Galactic extinction 
(jCarderri et al.lll989f) in the UV and NIR regions. 

For the idealized case considered in this work, a frac- 
tion of the photons scatter inside the dust shell and the 



r / r \ 
At r ~ - ~ 4 • — r- days. 



delay time induced by the multiple scattering will scale 
with r, 

10 16 cm7~~ J ~" (2) 

Since the typical rise-time of SNe la lightcurves is about 
3 weeks, and the fall-time is somewhat longer, we expect 
significant effects on the shapes of optical lightcurves 
for ?'J>10 16 cm for increasing optical depth, r(A). On 
the other hand, for time scales At r ^> 1 month, i.e., 
r S> 10 17 cm, we expect the time perturbations to de- 
correlate with the lightcurve shape parameters currently 
used, since these are designed to capture smaller time 
scale effects. 

6. LIGHTCURVE SHAPE 

Existing Type la supernova lightcurve fitters 
(IRiess et all [l99l IGoldhaberl 12001 IWang et al.l l200l 



Guvl|2005HGuv et al.ll2007t Uha et al.l[2007tlConlev eTaLl 
2008 : lBurns et al.ll2011| ) use different approaches to take 
variations of the optical lightcurve shape into account 
for standardizing SNe la. To estimate the effect on the 
lightcurve shape due to circumstcllar dust, we once again 
consider the two extreme scenarios: (a) = rwo ~ 
and (b) = 0.95 • r, introduced above. 
Fig. [2] shows the time-delay distribution in the B band 
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from Monte Carlo simulations as a function of the outer 
shell radius r for cases (a) and (b) where the LMC-type 
dust density was chosen to cause an observed color excess 
E(B-V) « 0.12. 

The far left panel of this figure is closely related to 
the left panel of Fig. Q] except that the dimensionless 
property D/r has been converted to days by assuming 
different physical shell radii. First, we note that for very 
large outer radii, photons that scatter on dust grains 
have a fairly flat delay probability with respect to the 
unscattered ones in the first (zero) bin. We also note 
how the effect that certain photon paths are less likely 
for a thin shell, as discussed above, propagates to a sharp 
feature in the delay time distribution. 

Also shown in Fig. [2] is how the delayed photons per- 
turb the pristine B band lightcurve shape. This reveals 
how the size of the outer radius, affects (primarily) the 
post-max shape of the lightcurve. Notably, the largest 
deviation from the pristine shape can be seen in the tail, 
more than a month after lightcurve maximum, typically 
beyond the time range investigated for most reddened 
supernovae. 

6.1. The fall-time 

The most straight-forward method of quantifying post- 
maximum lightcurve shape is in terms of the magni- 
tude difference between maximum and an arbitrary num- 
ber of rest-frame days past maximum in a given rest- 
frame passband. T his method of m easuring the fall-time 
was introduced by iPhillipsI (|1993t ) who chose the rest- 
frame B band and 15 days past maximum respectively, 
Ami5(B). The simulation results of the evolution of the 
Am\^{B) parameter for a few different CS-scenarios is 
presented in the far right panel of Fig. [2] and in more 
detail in Fig. El 

For 10 16 <^r<S5 • 10 17 cm, Ami^B) decreases with r. 
However, for larger radii, the late photons populate the 
tail of the lightcurve, i.e., and thus the original pristine 
lightcurve fall-time at earlier epochs is gradually recu- 
perated, while a plateau is built up at late times. The 
significance and the variation of the lightcurve tails in 
the examples shown suggests further studies of the fall- 
time, also at epochs beyond day 15. This is illustrated by 
the right panel of Fig.[3l where two different radii for the 
same dust configuration can give rise to the same value 
of Ami5(.B), but different values at 35 days past maxi- 
mum, Am35(i?). For very large radii, both the original 
values Aniis(B) and Am^(B) are recovered. For these 
extreme radii, the imprint from the CS dust will be found 
for even later epochs. 

Another interesting effect occurs for the dust shell sce- 
nario (b) as shown in Fig. [2] The fall-off of the r = 
10 18 cm curve is initially steeper than the r = 10 17 cm 
curve, but at day ~ 30 past maximum, this changes 
and the r = 10 18 cm curve flattens out and cross the 
r = 10 17 cm curve, while the latter continues to fall. 
This qualitative behavior can also be seen in real data: 
two SNe la normalized to the same maximum brightness 
show different relative fall-off for different epochs, as seen 
in Fig. |U 

6.2. The rise-time - fall-time relation 

The impact of the time-delay effect from scattering on 
CS material is expected to have a stronger impact on the 



fall-time than on the rise-time of the lightcurve. The rela- 
tion between the SNe la rise time and Arrii^B) has been 
studied both for high-qu ality data of well-sampled nearby 
SNe la (jStrovinkl I2007T ) . as well as in a statistical ap- 
proac h for a much large r SDSS -II data set (jHavden et al.1 
[2010^ . While [St T rmdhH (pOOl found double-peaked dis- 
tribution for the relation between ris e and fall time, this 
was not seen by (|Havden et al.|[20To() . Fig. [2] reveals that 
different CS-dust scenarios could give rise to a relation 
between rise and fall time, even when a single pristine 
light source is assumed. In Fig. [5] we show this relation 
for a few different s cenarios, together with the results 
from lStrovinkl (|2007l ). 

The absolute values of any lightcurve parameters de- 
termined from our simulations will depend on the cor- 
responding values of the pristine source. Since the 
lightcurve parameters of the H07 template will roughly 
reflect the mean values of the SN la population, our sim- 
ulations will only be able to produce lightcurves with 
longer rise times and smaller Am\^{B) than average. In 
the CS scenario it is reasonable to assume that a "naked" 
SN la template have shorter generic rise time and a faster 
fall-time, in order to compensate for this we have also 
applied an add-hoc shift to all curves in Fig. [5] in order 
to match the data points. Further, assuming that the 
generic B band shape of a "naked" SN la is similar to 
the H07 template that we use here, any modifications to 
the relative lightcurve properties for different dust sce- 
narios is likely to be of second order, and the general 
trend and size of the effect should still be valid. 

The result from Fig. [5] is that for SNe within a narrow 
range of observed color excess, CS-dust could naturally 
give rise to a bi-modal distribution in the relation be- 
tween rise and fall time for a small sample. However, if 
the accepted color range is extended, the gap between 
the two populations will be filled. More generally, it is 
intriguing that the rise-fall time scatter observed in nor- 
mal SNe la is compatible with our simulations. 

6.3. Time of maximum for different wavelengths 

Up to this point we have focused on the effect of cir- 
cumstcllar dust on the rest-frame B band, the wave- 
length range that is traditionally used for doing SN cos- 
mology. Since the scattering and absorption cross-section 
decreases with wavelength, we expect the delay distribu- 
tion to affect bluer pass-bands more. One way of study- 
ing this is to investigate the impact on the SN rise-time 
for different filters. However, this is related to comparing 
the time of maximum in different bands, which is also a 
property that can be measured quite accurately. Fig. [6] 
shows the effect on the difference of the time of maximum 
between the B and the / bands, t^ lax — £„ lax , for various 
CS scenarios. The general trend in the model predic- 
tion is that the time between the / band and B band 
lightcurve maxima increases for slow dcclincrs. Further- 
more, the effect is quite sensitive to optical depth, i.e., 
to the color excess. 

7. INTRINSIC COLOR VARIATIONS 

G08 explicitly describes both how circumstellar dust 
reddens the pristine source and derives this reddening 
law. However, as we have seen due to the time delay, 
mixing of photons originating from different epochs is 
expected. Since the delay time, At, depends on the shell 
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Figure 2. Simulated effects on the SN la B band lightcurve from multiple scattering by circumstcllar dust for three different radii: 
(blue, red, green) = (10 16 ,10 17 ,10 18 ) cm. The upper row shows results for r; = 0.0 while the lower is for r, = 0.95 ■ r. A reddening of 
E(B — V) ~ 0.12 has been assumed for both cases. The far left column shows the photon delay distribution in the B band. The middle left 
column show the corresponding B band light curves after the source (black) SED has been convolved by the delay distributions. The middle 
right column gives the B lightcurves normalized by source lightcurve at maximum. The far right column shows the B band lightcurve 
Amis (B) f° r different values of circum stellar shell radii. 



radius, the observed color for a fixed dust amount will 
have an intrinsic scatter with a magnitude depending on 
the distribution of the shell size, r. 

Fig.[7]shows how a few different CS-scenarios could af- 
fect both the lightcurve shape and the color at B band 
maximum. As with previous results we can see that the 
effects of both lightcurve shape and color increases with 
the amount of CS-dust. There is also degeneracy be- 
tween different CS scenarios, so no trivial correlation be- 
tween fall-time and color excess is expected in the CS- 
scenario. The CS-model thus naturally accommodates 
an intrinsic color scatter in moderately reddened super- 
novae, Qe( b-V) ~ 0-05 — 0.1 mag, in agreement with ob- 
servations (jNobili fc Goobarl 120081: iKessler etall I2010D . 



7.1. Time dependent color excess 

One of the salient features of the scenario in G08 is that 
blue photons scatter more than red photons. This implies 
that blue photons will, on average, be more delayed by 
multiple scattering since they have a longer path-length 
before leaving the circumstellar dust environment. This 
implies that the color excess due to the dust component 
will acquire time dependence. With respect to a time 
averaged color excess, the measured color will be redder 
pre-max turning bluer after max, as shown in Fig. [5J 
This effect should provide the most stri ngent bounds on 
the du st shell sizes allowed by the data. lYamanaka et al.l 
(|2009f ) show in their analysis of SN 2006X that once they 



correct for the color excess around maximum light, the 
V — R color is about 0.3 mag bluer than the normal 
unreddened SN 2003du 2 months after lightcurve peak, 
in qualitative agreement with the CS model prediction. 

8. PERTURBATIONS TO THE LIGHTCURVE 
SHAPE-BRIGHTNESS RELATION 

Since lightcurve shape is the primary parameter to 
standardize type la supernovae as distance indicators, 
it is important to investigate how the relation between 
lightcurve shape and brightness may be affected by the 
presence of circumstcllar dust. As noted in Section [31 
an intrinsically brighter explosion will deplete more of 
the surrounding dust, emphasizing the fainter- redder re- 
lation. For shell radii r>10 16 cm, the optical lightcurve 
shape will broaden the lightcurves potentially influenc- 
ing the empirically derived brighter-broader relation. Is 
it possible to obtain a correlation between brightness and 
lightcurc shape from the interaction with the circumstel- 
lar environment? Although possible, it invokes a some- 
what contrived scenario: For shell sizes r ~ 10 15 — 5 • 10 17 
cm, if the outwards velocity of the circumstellar mate- 
rial, i.e., the white dwarf wind velocity, correlates with 
the amount of 56 Ni powering the explosion, brighter SNc 
would have large r and become slow dcclincrs. Wc are 
not aware of any supernova models that have predictions 
in this respect. 

More generally, if the lightcurve shape is a combina- 
tion of effects, involving both physics of the progenitor 
system and the material around the white dwarf, then 
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Figure 3. Left: evolution of the lightcurve width expressed in terms of Amis(B) for cases (a) (dashed) and (b) (solid) for different 
amounts of dust. Right: the fall time in terms of Amis(B) and Am35(B) as a function of the shell radius. 
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Figure 4. B band lightcurves of 2006ax and 2006X from 
the Carnegie Supernova Project with Ami^(B) = 0.99 and 
Ami^(B) = 1.1 respectively. The solid lines are smoothed splines. 
While the fall-off between maximum and day 15 is steeper for 
2006X than for 2006ax, this is not true between maximum and 
for e.g. day 30. The same qualitative appearance can be seen for 
the CS model in scenario (b) in Fig. [2] if the shell radius is varied. 

one would expect a residual dependence of brightness 
on lightcurve shape, even after the main component has 
been calibrated out. In particular, we note that the opti- 
cal lightcurve shape beyond 30 days after B band peak, 
will be noticeably affected for a wide range of shell sizes. 

9. THE SECONDARY "BUMP" IN THE / BAND 
LIGHTCURVE 

The type la / band lightcurve exhibits a secondary 
maximum approximately 15-30 days after the B band 
lightcurve maximum which, in turn, typic ally happens 
withi n 2 days from the primary / band peak (|Nobili et al.l 
l2005f ). The time gap between the two I band peaks, as 
well as the relative strength of the secondary peak has 
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Figure 5. Rise-time vs fall-time imposed by different CS- 
scenarios. The colors represents different amounts of CS-dust and 
different points on the individual curves represent varying dust 
radii. The plot is for a spherical scenario, i.e. r; = 0. Different 
values on n give rise to similar size of the curves, but with slightly 
different shapes. 

been shown to corr elate with the B band lightcurve shape 
(jNobili et al.ll27)05ri . 

The implication of the imposed time-delay from scat- 
tering on circumstellar dust on the I band peaks is shown 
in Fig. IH1 We note that for a constant shell size we ex- 
pect the brightness difference (upper row) between the 
peaks, as well as the significance of the secondary peak, 
to decrease with an increased amount of dust. The effect 
comes from the fact that photons from the first peak are 
shifted to both fill the gap between two peaks and em- 
phasize the second bump, and is therefore also correlated 
with B band fall-time and E(B — V). 

If the shell size is allowed to vary, the usual degeneracy 
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Figure 6. The difference between the time of maximum in the B 
and / bands for different amounts of CS LMC dust. The curves 
are shown for a dust sphere, case (a). For the one of the curves, a 
few different outer radii have been marked for guidance. 
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Figure 7. The evolution of Amifj(B) and the color excess at B- 
maximum as a function of the shell radius for cases (a) (dashed) 
and (b) (solid) respectively. 

between shell size and the amount of dust is obtained for 
the observables. Despite this degeneracy, it is interesting 
to note that for a broad range of dust properties and shell 
sizes, a relation between the / band lightcurve properties 
and E(B — V) is expected. 

10. SPECTRAL FEATURES 

iFolatelli et all (|2010f ) found that different fitted red- 
dening corrections were obtained for their data set de- 
pending on whether the reddest SNe were included in the 
fit or not. For the full data set they obtained Ry = 1.7, 
while they fitted Ry = 3.2 if the high- reddened SNe were 
excluded. In contrast to this, lAmanullah et all ([2010D 



found that redder SNe in the Union2 sample prefer a 
higher color correction than bluer SNe and speculate that 
the color correction may very well be more complex than 
a simple li near relation . 

Further, Wang et al. (2009j) found a bi-modality in the 
color-magnitude relation over a broad color range, with 
the two groups preferring different values of Ry inde- 
pendently of the color. They also observed a correlation 
between the velocity of Si II(A6355) and the fitted value 
of Ry, concluding that objects with high velocity tend 
to prefer a lower value of Ry than SNe with low velocity. 

One possible explanation for the bi-modality could be 
that the color of the low-i?y SNe primarily originates 
from CS dust reddening, while the high-i?y objects are 
dominated by extinction from interstellar dust. 

The time-delay of photons from CS dust scattering is 
expected to blend the spectral energy distributions be- 
tween different epochs, i.e. a spectral feature for a given 
epoch will effectively be a superposition of the feature for 
all previous epochs. Since time-delay distribution varies 
with wavelength, the amount of blending will also de- 
crease with wavelength. Furthermore, since the minima 
of the supernova absorption features trace the velocity 
of the receding photosphere during the first weeks af- 
ter the explosion, we expect that delayed photons ar- 
riving around maximum light will originate from a re- 
gion further out in the photosphere, i.e., at higher ve- 
locities. However, our simulations indicate that these 
effects are small (see Fig. [T0|) . typically accounting for 
velocity variations of a few hundred km s , i.e., com- 
patible with the velocity scatter between measurements 
of normal SNe la, but not enough to explain significant 
outliers like SN 2006X (jYamanaka et al.ll2009fl . 

11. SUMMARY AND CONCLUSIONS 

Dust layers surrounding the progenitor systems of 
SNe la could help explain the observed reddening law. In 
this work, we have performed Monte Carlo simulations 
that indicate that the scenario with circumstellar dust 
would also perturb the supernova optical lightcurves in 
a manner that resembles empirical findings, e.g., an "in- 
trinsic" color variation of <Je(b-v) ~ 0.05— 0.1 arises nat- 
urally in our simulations. Since multiple scattering can 
only introduce a delay in time, the net effect is a broad- 
ening of the lightcurve. While the broadening increases 
monotonically with the density of dust, the relation with 
the dust shell size is less straight forward: there is an 
increase for outer radii about ~ 10 17 cm, while for even 
larger radii, the late photons populate the tail instead, 
and the lightcurve width is closer to the original one. 
We conclude that well sampled multi-band lightcurves of 
near-by SNe la reaching 2-3 months after lightcurve peak 
would be critical to probe the existence of dust in the 
circumstellar environment. Our simulations also suggest 
that for large optical depths, the lightcurve shape per- 
turbation should be notably different before and after 
maximum. Furthermore, unlike the case for reddening 
by interstellar dust, color excess should be epoch depen- 
dent: redder earlier on and bluer after maximum. 

A strong case for the CS scenario would be if the ob- 
served lightcurve perturbations in the optical could be 
combined with a detection of the emission at longer wave- 
lengths from the heated CS materia l. A m odel for this 
emission was constructed by iDwekl (|1983f ) for circum- 
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Figure 8. The imposed color for different epochs. The upper row shows the B — V color evolution as a function of time, for different 
amounts of dust (left) and different radii (right). The lower lower row show the color difference between 35 days after B band maximum 
and maximum brightness for different amounts of LMC dust. This is shown versus the B band fall-time (Ami5(B)) and the color at 
maximum in the left and right columns respectively. The simulated relation between the imposed lightcurve shape, The dust radii for a 
one of the CS realizations have been plotted along the curve for guidance. All results are shown for a CS scenario (a) with r; = 0.95 • r. 
The legend describes the color excess the different CS scenarios would give rise to if the time-delay effect on the color is ignored. 



stellar dust around Type II SNe. Although the exact 
properties of the emission is highly model dependent, it 
is generally expected to last after the optical light has 
faded. 
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APPENDIX 
A SPHERICAL DUST SHELL 

A simple model of a spherical dust-shell with radius, r, and inner radius, rj = R- r, where < R < 1 is illustrated 
in Figure IAT1 Also shown in the figure is a photon leaving the sphere (solid line). The current position of the photon 
is defined by the vector, f = (x, y, z), and the next position, f nex t, is given by f next = f Q + frj, where fn is the 
displacement vector fu = (dx,dy,dz). 

Photon traveling inside the shell — The impact parameter, i, is the minimum distance to the center of the sphere for 
any given photon-path, and can be defined as. 

. _ J 7*o • sin 7 if < 7 < | , 
l ~\r if f <7- 
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Figure 10. The the impact on the spectrum at maximum from CS dust. The upper panel shows the evolution for different shell radii, 
assuming LMC dust that would give E(B — V) ~ 0.25 in the limit when the radius goes to zero. The lower panel shows the spectrum for 
different amounts of dust when the radius is fixed to r = 10 17 cm. The inner radius is fixed to n = 0.95 ■ r for all cases. 



Here, 7 is given by 

To ' H3 = ^o r D cos(7r — 7) = x • dx + y • dy + z ■ dz 

x ■ dx + y ■ dy + z ■ dz 



cos 7 = — cos(7r — 7) 




Figure Al. Illustration of a scattered photon (thin solid line) propagating through dust shell. The figure shows the intersecting plane of 
the sphere that is defined by the two vectors f o and ?d ■ 

If i < 7"j, the photon will cross the inner radius of the shell and then re-enter the shell. This will in turn extend the 
mean free path of the photon by the amount 2q = 1\fr\ + i 2 , which is the length of the path within r;. 

Path length — For each distance between interactions, td is added to the total distance traveled by the photon before 
leaving the sphere. For a scattered photon leaving the sphere, the last distance, s, traveled inside the sphere is given 
by 

s 2 = r 2 + r„ - 2rr cos (ir — 9 — 7) , 

where cos 9 = sin 7 • ro/r. The distance traveled by a scattered photon should be compared to the distance traveled by 
a non-interacting p hoto n, which means that the distance, d, between the surface of the sphere and the plane (thick, 
solid line in Figure IA1[) perpendicular to frj at distance r from the center must be added to the total distance, where 
d is given by 

Q 

l = 2r ■ sin - 
2 

7T 9 tt-9 

(fi = TT 



2 2 2 

n ■ 7T — 9 „ . n 

d = l ■ cos <p = 2r ■ sin - cos = 2r • sin 



2 2 



The total distance traveled, D, can then be summarized as 



jj _ ( r for non- interacting photons, (Al) 

~ \ J2 r D n + s + d for scattered photons . ^ ' 
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